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Abstract 

A detailed analytical and numerical analysis for the dislocation 
cloud surrounding a disclination is presented. The analytical results 
show that the combined system behaves as a single disclination with an 
effective fractional charge which can be computed from the properties 
of the grain boundaries forming the dislocation cloud. Expressions are 
also given when the crystal is subjected to an external two-dimensional 
pressure. The analytical results are generalized to a scaling form for the 
energy which up to core energies is given by the Young modulus of the 
crystal times a universal function. The accuracy of the universality 
hypothesis is numerically checked to high accuracy. The numerical 
approach, based on a generalization from previous work by S. Seung 
and D.R. Nelson (Phys. Rev A 38:1005 (1988)), is interesting on its 
own and allows to compute the energy for an arbitrary distribution 
of defects, on an arbitrary geometry with an arbitrary elastic energy 
with very minor additional computational effort. Some implications 
for recent experimental, computational and theoretical work are also 
discussed. 
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1 Introduction 



Topological defects, mainly disclinations and dislocations, play a crucial 
role in the understanding of the physics of many two-dimensional (2D) sys- 
tems pQ. Although computational techniques such as analytical methods, 
Ewald summation techniques [2], Monte Carlo simulations etc.. can be 
used in many different contexts, they become quite inefficient in other sit- 
uations, particularly with crystals having a boundary or lying on a curved 
background. In this paper, a new approach based on previous work by Se- 
ung and Nelson [I] will be discussed in detail. The computational methods 
developed allow to obtain the lowest energy configuration for an arbitrary 
distribution of defects in a crystal lying on an arbitrary geometry with an 
arbitrary elastic energy at zero temperature. The method is general enough 
to include Bravais lattices other than the triangular case. A discussion of 
the screening cloud of dislocations surrounding a disclination will be pre- 
sented. This problem will be used as a testing ground where the results may 
be compared with the analytical predictions from elasticity theory derived 
in this paper. 

An interesting experimental example where the dislocation cloud sur- 
rounding a disclination appear is given by colloidal particles crystallizing on 
the surface of a sphere (colloidosomes) [5]. As a consequence of the Euler 
theorem (5], an sphere must have a total disclination charge of 12. If the 
total number of particles forming the sphere is large enough, the ground 
state contains more defects than the 12 necessary to satisfy the Euler the- 
orem. Those additional defects are dislocations surrounding a disclination, 
as illustrated in fig. ^ arranging in the form of grain boundaries. Those 
grain boundaries have two very distinctive features 1) terminate inside the 
medium and 2) have a total disclination charge of +1. These features are 
intimately related to the geometry of the problem and should appear when- 
ever the Gaussian curvature is large enough 03 E]. Similar structures 
have been also observed in simulations of the Thomson problem (101 1111 112j . 
Clouds of dislocations (with the minus disclinations playing the role of plus 
charges) should also appear in crystals on a negatively curved background 
|13j . Those situations are relevant as two dimensional analogs of the frus- 
tration associated with the three dimensional tetrahedral packing ^1] and 
other surface physics problems ^3 • 

The problem of grain boundaries radiating out of disclinations is also im- 
portant for understanding certain aspects of the problem of two-dimensional 
melting. The KTHNY |171ll8llTU] scenario predicts a two stage melting from 
a crystal into an isotropic phase via an hexatic phase. This picture has been 
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confirmed by a large number of examples (see |2U1 |3j for reviews). Very 
precise numerical simulations [2J have found very good agreement with 
the predictions of the KTHNY scenario. Different inherent structures (IS) 
characterize each phase, and disclinations appear as forming grain bound- 
aries with very similar characteristics as the ones found in the Colloidosome 
problem. Similar structures have been found in recent simulations, and the 
defects have been characterized by 1/ f noise [221 • Recent experiments on 
plasmas |23] and colloids |241 125j . have found also a similar situation, with 
correlation functions in good agreement with the KTHNY scenario and grain 
boundaries with similar characteristics. Alternative melting scenarios have 
exploited the short-range nature of the stresses produced by grain bound- 
aries of dislocations [2f)j . 

The physics of Langmuir monolayers has received a lot of attention in 
recent years [23I2E]) but many important questions have not yet been clar- 
ified. It seems very clear that topological defects play a crucial role in 
problems such as melting or the collapse of monolayers [2H]. More sophisti- 
cated systems, such as sphingomyelin, where additional hydrogen bonding 
may be formed, show a much richer phase diagram 30 . Biphasic surfactant 
monolayers show additional defect structures (mesas) |31j . 

Radial grain boundaries radiating from a central disclination are also 
found in hexatic I* and crystal J* tilted liquid crystal phases |32j . The tilt 
is used to force a disclination and a pattern of radial grain boundaries (with 
5 arms) is observed. 

The methods described in this paper are also relevant for a restricted 
type of quantum dots |33j in which the density of electrons is is small and 
the external disorder potential is weak enough so that the electrons in the 
dot form a Wigner crystal, so called "Wigner Crystal Islands" [33j. The 
methods presented in this paper are relevant to this case as well, as only a 
change of the boundary conditions for the problem is required. 

Dislocation clouds also appear in many other problems such as, par- 
tially polymerized membranes [35[ I36| . Wigner crystals |3Z1 and in carbon 
nanotubes, particularly in the so called "Onion" rings, which are successive 
spherical layers of graphite [38] . 

The organization of the paper is as follows. A review of continuum re- 
sults is given in Sect. [2j The discretized approach is introduced in sect. 01 
and compared against the continuum results for isolated defects. The scal- 
ing relations satisfied by the energy are derived in sect. El Numerical results 
are presented in sect. Eland their universality with respect to the Lame co- 
efficients is discussed in sect. El Sect is a brief overview of the effect of 
curvature. We wrap up with some conclusions in sect. El Several technical- 
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Figure 1: Ground state configurations for a large and small colloidosome. 
For a large colloidosome, finite length grain boundaries radiate out of the 
disclination (from 5 ). 



ities are relegated to the appendices. 



2 Continuum Results 

The elastic energy of a continuum is given by [SHI 

F=\j ' d 2 r{2^ul p + \{u aa ) 2 ) (1) 
where the strain tensor is defined as 

a ^ 2 dxp dx a dxp dx a 

The quadratic term in the strain tensor is usually dropped, as it usually 
amounts to higher order negligible corrections to the total energy. Those 
terms cannot be entirely neglected if disclinations are involved [2] (For con- 
sistency, other energy terms quadratic in the strain tensor should also be 
included, but for the sake of simplicity this point will be ignored). The 
Young modulus and Poisson ratio are 

H + A 



K = A(i 



2/x + A 
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Using known analytical techniques, the strains for a dislocation and a 
disclination may be solved in linear order |39l For a dislocation with 
Burgers vector b = be x the result is jHU] 

Ux = ^ + |l Sin(20)) (4) 
The strains for disclinations of charge s are [1] 

u<t> = ^r<f> , (5) 

There is an arbitrary constant A s which will be determined in the appendix 
sect. \M 

Plugging Eq. 0] and Eq. [5] into Eq. ^ The energies for a dislocation of 
Burgers vector b is 

|b| 2 

-K (]n(R/a) + const). (6) 



8vr 

As it is shown in the appendix sect.EJ the constant A s is related to the two- 
dimensional pressure II. Following the same steps as for a dislocation, the 
energy for an isolated disclination of charge s = ^qi under pressure becomes 

F = —K R 2 + ^-R 2 . (7) 
32vr u 823 y J 

The energy dependence, growing as i? 2 , is energetically very costly. 
Grain boundaries of dislocations may screen out the disclinations, thus re- 
ducing the huge energy cost of an isolated disclination, if the angle of the 
grain boundary exactly compensates the missing or additional wedge (a 
multiple of ? for triangular lattice) caused by the disclination. From the 
geometric argument in fig. [2 the angle of grain boundary 9 is given by 

2i?sin(6»/2) = nb, (8) 

where n is the total number of dislocations all having Burgers vector b. 
Further assuming a constant spacing of dislocations within the grain D, the 
relation B = Rb/D holds and one gets 

2sin(0/2) = b/D. (9) 
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Figure 2: Relation between the aperture angle 9 of a grain boundary of 
dislocations and its total Burgess vector nb and radius of the crystal R. 



upon identifying 6 with the missing or additional wedge removed to form 
the disclination j^j, Eq. ED becomes 

2sm(-^) = b/D, (10) 

0771 

leading to an equation for the spacing D as a function of the number of 
arms m. A detailed analytical proof for this result using linear elasticity 
theory is provided in appendix El where the total energy of the system of 
a disclination and an m-arm grain boundary (see fig. [5J) is given by Eq. 1531 
There is also a linear term in R arising from two different contributions, the 
stresses of grain boundaries of dislocations, discussed in appendix^ and the 
core energies of the defects. The total contribution to the energy is then, 

(s - mi) 2 9 / , a , \ Kn a 2 , . 

F = 327r~ ^° + {"^ + H m lZ D R ' < U > 

where the function H is defined in the appendix Eq. |S7| The pressure II has 
been set to zero. The total number of dislocations is n = rn^. The core 
energy has been parametrized as E core = Kqo?c, where c is a dimensionless 
coefficient (but not independent of the elastic constants). 
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If the spacing D is given by Eq. I1U1 the leading R 2 term in Eq. ^2 is 
canceled and only the linear term in R , which we denote as /, survives 



f=(H(2^) + 4Trc^m^R, (12) 
For future reference, we quote the result for large m 

{= (r, (1 - k ^ + c ) mK °" 2 i- (13) 

At perfect screening m|j = <2i§ and the energy grows logarithmically with 
m. All previous results assume an infinite system (large R). How large R 
must be in order that the infinite radius result hold, will be estimated next. 

The interaction energy of two grain boundaries decays exponentially fast 
as a function of their mutual distance [10], with a decay length 

V = (14) 

D being the distance between dislocations within the grain (fig. EJ). 

dislocations that are a distance L from the disclination will be invisible 
to the dislocations in other arms if 



2Lsin(vr/m) >> 2A 9 . (15) 

therefore from Eq. [21 and Eq. EH it follows 

3 

L » -m 2 b. (16) 

7T 

For R > L the infinite result Eq. ^Jwill hold. The dependence of Eq. H6l on 
m 2 points out that the behavior Eq. ^]for large number of arms will not be 
reached until the radius R is very large. 



3 Discretization of the Elastic Free Energy 

We consider a 2d triangular lattice of monomers with lattice constant a. 
The actual position of the (n, m) monomer is described by Yt n ^ m \ and it may 
be decomposed as 

r(n,m) = nei + me 2 + u^ m) , (17) 
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where is the strain at point b and ej, i = 1,2 define a basis of the Bravais 
lattice. For a triangular lattice it is 



ei = e x (18) 

1 ^3 
-c r H 

2 2 

The discrete elastic free energy that we will use is 



e 2 = 7: e a; + — e y ■ 



F = \Y. (Ir 6 d - I) 2 + - £(5 " i^m^ 1 . ) 2 • (19) 

<6c> 6,c ^ l r bc||r b c+l| 

The summation in the first term runs over links defined by vertices a, b, and 
in the second term the sum runs over all nearest neighbors of point b. Other 
discretizations are also possible but do not modify the results as it will be 
discussed. 

In order to proof the equivalence of Eq. ^3 with the continuum result 
Eq. one must define the discrete derivatives first. The derivatives of a 
general function /(r) are only defined along the discrete direction of the 
triangular lattice, so a general derivative is precisely defined from 

d x f(n) = /fo + ^i) -/('») (20 ) 
a 

o f( , fin + ae 2 ) - f(r b ) - \{f{v b + oei) - f(r b )) 
yHrb) ~ aV3/2 

Using the previous definition, the discrete metric becomes 

gij(b) = diT b djY b . (21) 

The distance between two nearest neighbor points is 

|r a -r b | = (l + 2^-4 6 4 6 ) 1 /2 ) (22) 

where the discrete strain tensor is defined from the derivatives Eq. |2U] by 

Q/3 2 dxp dx a dxp dx a 

Expanding Eq. [221 using Eq. ED derived in the appendix, the first term of 
the discrete energy is identical to Eq. ^ with 

X ^ 
= ~ e 

M = ^e. (24) 
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The same expansion to the second term using the relations Eq. 1601 1611 1621 
from the appendix lead to 



A = 9^ ff 
4 

/i = -9^6. (25) 
The elastic constants of the two terms combined are then 

A = ^(e + 9a) 

H = ^(e-9a), (26) 

and therefore, Eq. ED provides a suitable discretization of Eq. ^ with arbi- 
trary elastic constants. It should be recalled that higher order terms in the 
displacement are dropped to lead to Eq. |2T)1 

In this paper, the particular geometry that will be used consists of a 
plus or minus disclination (a pentagon or an heptagon respectively) at the 
center of the crystal. The total number of monomers forming the crystal is a 
function of the linear size R and depends on the total number of dislocations. 
If only a center defect of charge qi is present, the total number of monomers 
is 

M = (6 — qi)—^— + 1 , (27) 

Even with additional dislocations, the total number of points will still grow 
like (6- qi )R 2 . 



3.1 Energies of single defects 

The energies of isolated dislocations and disclinations will be compared 
against the analytical predictions derived previously. This will provide a 
benchmark to the present approach. 

In fig. |21 the minimum energy results of a configuration containing a 
dislocation for different values of the Lame coefficients. A fit to the form of 
Eq. El yields 

-4-2 = 0.0395(2) log (r) (28) 

in agreement within of 1% with the analytical result = 0.03979. The 
small deviation may be attributed to the neglect of higher order terms in the 
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Figure 3: Results for the energy of an isolated dislocation as function of R 
for different values of the elastic constants. 
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Table 1: Core energies of dislocations as a function of Lame coefficients 



continuum calculation. Core energies may be computed from the intercept 
in fig |21 The results are summarized in tabled It should be noted that core 
energies for crystals with the same Young modulus are different. 

The results of the same analysis for single disclinations are shown in 
table |2j The most accurate determination yields 

^=0.00785(1), (29) 

both for plus and minus disclinations (Results for A = /i were first obtained 
with less accuracy by ; 4] ) • Although the coefficient is off from the analytical 
result Eq.[7|by a significant amount, it is remarkably universal as a function 
of the elastic constants. The energies are also the same for both plus and 
minus disclinations within a 0.1% accuracy (more accurate results do show 
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Figure 4: Plot of function E{r^) describing the relative difference of the 
analytical solution and the numerical result. 
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Table 2: Coefficient (Eq. l29|) as a function of the Lame coefficients (/i, A) 

that sevenfold defects have a marginally lower energy). There is, however, 
a term that grows linearly with R in the discretized energy. This term is 
quoted also in table |2j 

It is also instructive to compare the continuum strain solution r A Eq. |S] 
with the configuration from the discrete calculation. The difference will be 
quantified from the function 



E(r b ) 



N + \ t a\ 



(30) 



and it is plotted in fig. El The relative error is always below 0.01%, becoming 
as low as 10 _6 %. 

4 Scaling Relations for Grain Boundaries 

The degrees of freedom defining a grain boundary of dislocations are (see 
fig.©: 



• m: number of arms. 
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Figure 5: Degrees of freedom of a grain boundary. 



• Igb'- Length of the grain boundary from the central disclination. 

• angle of the grain with some specified crystallographic axis. 

• D: spacing of dislocations within each grain. 

The spacing D does not generally need to be restricted to be constant. The 
linear size of the system is R. 

It is convenient to define the dimensionless variable 



which by definition satisfies < x < 1. x = 1 implies that grains reach the 
boundary while x = implies no dislocations. Eq. will be generalized 
(for zero temperature) by assuming that for intermediate values of x the 
following scaling behavior holds 



where as x — > Q approaches Eq. [Tj This scaling law however, may break 
down if the sub-leading linear terms in R become comparable. The condition 
expressing this situation is given by (see Eg. I12JI . 



X = lgb/R 



(31) 



E = K R 2 Q(x,m,D,V) 



(32) 




(33) 
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If the term square in R vanishes, this equation must hold for some x c and 
scaling will break down for all x > x c . 

If the square term in R does not cancel, then the scaling relation must 
hold for large R as well. That will be the case for spacings D not satisfying 
Eq. I1UI The Q function should exhibit a minima roughly at the critical 
x m where the additional angle added by the dislocations compensates the 
missing or additional angle by the disclinations. The critical x m is given by 



with the additional constraint x m < 1. Therefore the Q should exhibit a 
minima for x m . 

The angle ^, defined as the angle of the grain with respect a crystallo- 
graphic axis, for an m-grain boundary commensurate with the p-fold sym- 
metry of the central disclination will be constrained to 



A straight-forward analytical calculation shows that Q should be indepen- 
dent of It will be shown that this result ignores the constraints imposed 
on the Burgers vector by the underlying lattice. 

5 Numerical study of grain Boundaries 
5.1 Some computational details 

The calculations have been done by relaxing an initial configuration con- 
sistent with the given distribution of defects using the conjugate gradient 
method. The numerical accuracy was tested by checking the convergence of 
the final results as a function of the tolerance error in the algorithm. When- 
ever different initial configurations consistent with the given distributions of 
defects were tried, the final result was found to be identical. The energies 
were computed to a seven digit precision or more. 

For each value of the parameters, results were obtained for linear sizes 
ranging from R = 10 to R = 200 corresponding to typical volume sizes 
from 247 to 140000 monomers. Those lattices are, in many cases, larger 
than some experimental systems available. Free boundary conditions were 
incorporated by allowing the system to reach its natural extend without any 
external constraint. 

The code has been written in C++ using objected oriented design. This 
provides the flexibility of incorporating any potentially new geometry, dis- 
tribution of defects or discretization energy at any future time with a very 




(34) 



< * < — . 

m 



(35) 
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minute effort. The code has some limitations too, one of them being that 
non-integer spacings (cases were dislocations within the grain should be 
slightly non-constant to mimic fractions of the lattice spacing) have not been 
implemented effectively. Typical minimizations of a system with R = 100 
(total number of monomers M = 25000) take 9 minutes at a 10 -10 precision 
and 5 minutes at 10 -5 precision in a Dell 1.80 GHz dual Xeon procesor 
running Linux Red Hat. Further computational details will be presented 
elsewhere. 

5.2 Scaling as a function of l g b and D 

A plus disclination will be placed at the center of the crystal and grain 
boundaries will be fixed to have angle Vl/ = 7r/5 and number of arms m = 5. 
The energy will be investigated as a function of the parameters l g b and D. 
The elastic constants will be chosen as A = [i. 

The energy of relaxed configurations is plotted in fig. El as a function of 
the scaling variables, the Q function is plotted for a plus-disclination and a 
m = 5 grain boundary with spacing D = 5. The plots follow very well the 
ansatz in Eg. 1321 with obvious deviations for values of x closer to 1. 

In order to determine the coefficient of the linear term in R, defined as 
/ in Eq. 1121 it will be assumed that the points on fig. El for x ~ 1, which 
clearly do not scale, are described by Eq.^J The results are plotted in fig.0 
and one obtains for the /-coefficient 

/ = 0.081(5) . (36) 

Since we are not looking for a high precision value for the /-coefficient, two 
values for l g b were included for any given R. This provides strong evidence 
of the robustness of the fit. The theoretical value of the / coefficient arises 
from a grain boundary contribution Eq. [SHI and a core energy which may be 
estimated from a single dislocation (given in table^). The two contributions 
give 

/ = 0.063(Grain Boundary) + 0.029(Core Energy) = 0.092 , (37) 

considering the approximations involved (basically linear elasticity theory 
and non-interacting grains), the agreement with Eq. 1311 is acceptable. 

The value obtained for / may be cross-checked by assuming that the 
scaling behavior will break down when the value for the linear term is, say 
one third, of the square term. The intersection of the straight lines in fig. El 
provide a visual solution to Eq. ESI and an estimate of the critical Xq clS EL 
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Figure 6: Plot of the scaling function Eq.|22for D = 5, m = 5, Vl/ = vr/5 and 
sizes i? = 10 — 170. The dashed line is a fit with an extrapolation to x— > 1. 
The straight lines provide a visual solution to Eq. [2S1 an d an estimate of the 
critical function of R. 
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Figure 7: Plot of the fit Eq. HUto the results for R ~ l gb . 
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function of R in reasonable self-consistent agreement, particularly for large 
values of R. 

Three regions may be clearly identified from fig. |HJ 

• (x ~ 0): The energy is dominated by the strains of the central discli- 
nation and the effect of adding additional defects is negligible. 

• (0 << x < x c ): The inclusion of more defects dramatically lowers the 
energy. 

• (x > x c ): The core energy contribution sets in and scaling breaks 
down, the energy grows linearly with R, as apparent from fig. 

It should be noted from fig.EJthat if the intermediate region could be extrap- 
olated to to x ~ 1 before the core energy terms would become noticeable, 
the energy would go to zero as (x — l) 2 , that is, independent of the linear 
size R. 

A typical relaxed final configuration is shown in fig. 1121 The plot is only 
for a small region around the central disclination, but it clearly shows that 
for regions away of the defects the triangles forming the crystal are almost 
equilateral (unstrained) . 



5.3 Scaling function for non-optimal D values at fixed m and 

The Q function for D values significantly different than D cr n = |vr are 
plotted in fig. (D = 3) and fig. (D = 7) . 

For D < D cr n the scaling function has a minima as a function of x 
in reasonable agreement with the value predicted by Eq. 1341 As it is very 
apparent from fig. |5] the energy grows quadratically (compare with fig. |5J) 
with the system size, even in the presence of the grain boundaries. The fit 
gives 

= 0.0055(4) . (38) 

This result should be compared with the theoretical estimate Eq. [221 for 
n = 0, D = 3, m = 5 

Although the result is not too different, the quantitative agreement is not 
very good. This may be attributed to assuming linear elasticity, which as 
already seen is not very accurate if disclinations are involved. 
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Figure 8: Plot of the scaling function Eq. EUfor D = 3, A = [i, m = 5 and 
* = f for sizes R = 10 - 200 (A = fi). 
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Figure 9: fit Eq. [TT]to the values R~ l gb . 
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The energy for D = 7 does not exhibit a minima because that should 
appear for values of x > 1. The total energy for D = 7 is significantly larger 
than for D = 5, as it becomes apparent from the relative scale of the y-axis 
of the plot. The theoretical estimate Eq. [231 for II = 0, D = 7, m = 5is 

0.0011 . (40) 



KnR 2 



The result in fig. 1111 approaches this limit, although the statistics are not as 
good as in for D = 3. 



5.4 Dependence on the orientation of the grain ^ 

The dependence on the angle ^ is plotted in fig. ^] for angles Vl/ = 0, jq, |. 
Typical final relaxed configurations are shown in fig. ^] = 7r/5) and 
fig- EH = 0). The final relaxed configurations for ^ = ^ are very regular 
while the dislocations for the angle $ = display a rather jagged pattern. 
Provided that the total Burgers vector is zero, the only difference in the 
energies arises from the grain boundary terms. Therefore, the energy differ- 
ence in the results fig. EH should be attributed to the constraints induced by 
the lattice to the Burgers vector. This point might be substantiated numer- 
ically with a more comprehensive calculation, but this has not been done, 
partly because the angle is not well defined for values of m different than 
5. The previous considerations reflect that the optimal angle for m = 5 
grain boundaries is given by = n/5. 



5.5 Dependence on the number of arms m 

The analysis for different value of m will be restricted to the optimal spac- 
ing values Eq. |§] and to the more relevant case x ~ 1 (l g b ~ R). Some 
representative relaxed configurations are shown in fig. I14H17I Similarly as 
it was found in the investigation of the dependence on \& subsect. 15.41 the 
constraint that dislocations can only be oriented along directions defined by 
the triangular lattice is the origin of additional frustration, leading to jagged 
arrangements which follow the ideal orientations only approximately. 

The results corresponding to the energy are shown in fig. 1181 For finite 
radius, the smaller values for m are clearly favored. As R becomes large the 
values for the energy become degenerate (within the numerical accuracy) 
in m for m within the range m = 2 — 5. An evaluation of the H function 
in Eq. El points out that the m = 2 result should have the lowest energy. 
However, the difference in energy for arms m = 2 — 6 are small compared 
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Figure 13: Final configurations for m = 5, D = 5 and = and A = fi. 
Results correspond to R = 90. 
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Figure 14: Relaxed configurations for m = 3, D = 3 and A = fi. Results 
correspond to R = 90. 




Figure 15: Relaxed configurations for m = 4, D = 4 and A = /i. Results 
correspond to R = 90. 
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Figure 16: Relaxed configurations for m = 7, D = 7 and A = \i. Results 
correspond to R = 90. 




Figure 17: Relaxed configurations for m = 10, D = 11 and A = fi. Results 
correspond to R = 90. 
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Figure 18: Plot of the energies for l g f, ~ R as function of R for A = \i. m = 3 
(square) m = 4 (cross) m = 5 (circle) m = 7 diamond and m = 10 (triangle 
down) . 



with the core energy contribution, and at this precision, contributions from 
second order elasticity theory which have not been included may become 
important. 

For larger m the results show a larger energy, which is in qualitative 
agreement with the logarithmic dependence given by Eq. 1131 Numerical 
results for large m become increasingly difficult because it is difficult to direct 
the dislocations along the correct directions and also, because of the small 
number of dislocations per arm involved. A study including larger volumes 
than the ones performed here is necessary for more rigorous results. Finite 
size effects have been predicted to be negligible quadratically as a function of 
m (see Eq.ll2|). The convergence is roughly consistent with that dependence. 

5.6 Scaling collapse for minus disclinations 

Results for minus disclinations have not been computed with the same ac- 
curacy as for positive ones. Nevertheless, the same trends as for fivefold 
defects are observed, as shown in fig. 1191 The data collapses well to the 
assumed form Eq. I32( and the overall energy in the intermediate region is 
smaller than for fivefold defects. 
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Figure 19: Plot of the scaling function Eq. 021 for D = 7, A 



fi, m 



7 and 



for sizes R = 10 — 200. The dashed line is the universal function for 



a plus disclination. 
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Figure 20: Plot of the energy with the core energy contribution subtracted 
for l g b ~ R as function of R for m = 5 and different values of the Poisson 
ratio v. 
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Figure 21: Crystal on an spherical cap. No additional dislocations are 
needed. 

6 Universality of the Results 

The main assumption in the scaling forms for the energy Eq. |321 is that the 
only dependence on the elastic constants arises from the Young modulus. 
This is also true for all sub-leading terms, except for the core energy coeffi- 
cient c as it is obvious from tabled Therefore, we will present the results 
by subtracting out the core energy contribution, and the results should then 
become universal. 

Energy values for x ~ 1 plotted this way are shown in fig. l2Ul Except for 
very small systems, the universality hypothesis holds very well. Although 
cannot be conveyed from the plot, the energies for different Poisson ratio 
collapse to the same universal result with an accuracy less than 0.1%. This 
proves that the assumed form for the energy dependence is indeed universal, 
independent of the microscopic details of the lattice. 

This universality of the Q function is also remarkable since the present 
calculation goes beyond linear elasticity, not only through the quadratic term 
in displacements kept in the strain tensor Eq. I23[ but also through the higher 
powers of the strain tensor implicitly neglected in going from the discrete 
energy Eq. to the continuum result Eq. ^ Those non-linear elasticity 
terms have shown to be quite important when disclinations are involved, 
and yet, do not seem to disturb in any way the scaling form assumed for the 
Q function. It seems reasonable then to assume that the results presented 
are not only universal with respect the two elastic constants, but also with 
respect to higher order terms in the elastic energy. 



25 



7 The Spherical Cap 



Although the results presented in this paper are interesting on their own, 
much of the motivation for carrying out such a project arises as an effort 
to provide efficient computational tools for investigating crystals on frozen 
geometries. We will therefore present a brief outline of some preliminary 
results for that problem. This will also provide a broader perspective of the 
versatility of the method. 

Screening of disclination by grain boundaries is not the only mechanism 
available if disclinations are allowed to buckle out of the plane: Gaussian 
curvature without the need of additional grain boundaries is also a viable 
mechanism [IJ|H]. 

An spherical cap of aperture angle 7 will have a Gaussian curvature 

2 

K = ^7. Therefore, for large 7 the Gaussian curvature will suffice to screen 
out the disclination and no additional dislocations will be needed. The 
limit of vanishingly small 7 has been discussed at length in this paper. 
Disclinations are screened out by grain boundaries going all the way to 
the boundary of the spherical cap. Therefore, for intermediate values of 7, 
structures of grain boundaries interpolating within these two cases should 
be observed 

The result of a minimization for large 7 is illustrated in fig. 1211 It is found 
that, opposite to what happens for the flat case, additional dislocations 
actually increase the elastic energy of the system. One should notice that 
the triangles next to the boundaries are equilateral, which implies that the 
strains are very small (and so is the energy). A more detailed presentation 
of the results for crystals on spherical caps (positive curvature) including 
an investigation of the intermediate regime as well as a similar analysis for 
some minimal surfaces (negative curvature) will be presented elsewhere. 

8 Conclusions 

8.1 summary of the paper 

The advantages of the computational method presented are: 

1. Crystals with boundaries can be treated very efficiently. 

2. There are no long range interactions. An entire sweep over the whole 
system may be performed with a time proportional to the total volume 
of the system. 
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3. The calculations may be extended to additional geometries with a 
negligible additional effort (only introducing the coordinates defining 
the geometry). 

4. The convergence of the results (at zero temperature) is fast an stable. 

5. The results are very universal, valid for a wide range of potentials. 
Microscopic details only enter via the elastic constants. 

We provided a very detailed analysis for the problem of the dislocation 
cloud screening a disclination. It has been found that the system composed 
from a disclination charge s and m-radial grain boundaries of dislocations 
separated a distance D behaves as a single disclination with an effective 
charge 

b . . 

s eff = s- m— . (41) 

If s e j f = then the total energy of the system grows linearly with the system 
size, similarly as for an infinitely long linear grain boundaries composed 
of dislocations only. The systems exhibits a remarkable universality with 
respect the elastic constants (up to core energy terms) and higher order 
elastic terms. The analytical expressions derived have been compared with 
the numerical results, and when disclinations are involved the discrepancy 
occasionally may be as large as 25%. Linear elasticity theory does not 
provide accurate quantitative estimates for the energetics of some problems 
involving disclinations. 

8.2 Implications for Other problems 

We will briefly discuss some of the implications the results found in this pa- 
per have for the problems presented in the introduction. It has been shown 
that for very small systems, grain boundaries with the smaller number of 
arms have the lowest energy. This is a general result in agreement with other 
continuum calculations [Jj in the context of the sphere. The experimental 
results in Colloidosomes do show the same trends [3]. It would be very inter- 
esting to image larger Colloidosomes and analyze how the defect structures 
are changed, although that seems difficult in view of the equilibration times 
involved. 

It has also been shown that if disclinations appear in a crystal phase, 
grain boundaries of dislocations should follow. Isolated disclinations are al- 
most forbidden at zero temperature. At finite temperatures some properties 
of the grain will change, but for distances close enough to a disclination, 
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the huge strains that dominate the energy will make additional disloca- 
tions inevitable. Strings (m = 2 grains) or higher m grains should sur- 
round the disclinations, even at finite temperature. Most of these grains 
will have a non-zero disclination charge (positive or negative). This qualita- 
tive picture seems in agreement with recent numerical |21) an experimental 
results [241 123j , where the correlation functions do show agreement with the 
KTHNY scenario (and therefore with the unbinding of disclinations) while 
the typical snapshots of configurations in equilibrium contain grain bound- 
aries with non-zero disclination charge. Theories based on grain boundaries 
|26| ignore the possibility of having non-zero disclination charge from the 
very beginning. It has been shown in this paper that grain boundaries with 
non-zero disclination charge have energies comparable as grain boundaries 
of pure dislocations (with total disclination charge zero), provided the spac- 
ing within dislocations is fine tuned appropriately. A complete discussion 
of the temperature effects is obviously of great interest but it is beyond the 
scope of this paper. 

8.3 Outlook 

Some issues in this work need either further understanding or just higher 
precision data: 

• Dynamic defect distribution: A global minimization for dislocation 
positions and orientations would directly provide the minimum energy 
configurations. This is a difficult task since it has been shown that the 
energy has many almost degenerate local minima. 

• Non-integer spacings: Dislocations with non-integer spacings like D = 

= 5.23, which imply that every third dislocation must be separated 
an additional lattice constant, have not been implemented. This may 
be of some importance in connection with Eq. I1UI 

• The degeneracy of the energy as a function of m values should be 
further refined, since the accuracy of the results do not allow to dis- 
criminate for small values of m. 

• The effects of an applied pressure has not been investigated numeri- 
cally. 

The computational method presented can be used to compute the energy 
of any distribution of defects in a two dimensional crystal. Results are in 
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progress to investigate the effects of curvature and the interactions between 
disclinations. 

We hope that the considerable detail presented in this paper will not 
obscure the main results obtained, but on the contrary, will provide con- 
vincing arguments for the utility of the present approach. The code used in 
this paper will be made publicly available. It may be also requested from 
the author. 
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A Energy of a disclination 

The stresses of a single disclination at the origin may be computed from the 
Airy function of a disclination [1] 



Gs 
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(AR 2 + r 2 log(r)) 



(42) 



where s = ^qi,G = and A is an undetermined constant. 
The stresses for a disclination are given by 
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The stress along the radial direction is 



°rr = ^(A+i + log(r)) 



(43) 



(44) 



therefore, writing A = — 1/2+ — log(i?), the stress at the boundary R is 
o rr = IT, and IT is interpreted as the two dimensional pressure the material. 
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Plugging Eq. 03] m t° Eq. ^ the energy becomes 

s 2 o irll 2 9 

F = K Q R 2 + — -R 2 , 

32vr u 8B ' 



(45) 



where B = fi + A is the two dimensional bulk modulus. The result that the 
stress tensor Eq. 03] has an energy given by the expression Eq. 03] will be 
used quite frequently. 



B Analytical proof of the optimal dislocation spac- 
ing 

We will first compute the stresses generated by a finite grain boundary of 
length l g b, with dislocations within the grain separated by distance D. 

The stresses of a dislocation of Burgers vector b located at the y axis at 
point (0, z) are given by 
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(46) 
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where G = where Kq is the Young modulus. The total stresses will be 
computed from the methods used in dislocation pile-ups One obtains 
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A similar calculation yields 
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If the grain is now rotated and angle the stresses at distances r « l g b 
will be (where r = y/(x 2 + y 2 ): 
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Adding together m arms separated an angle — , one obtains 
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where the identity YaLi «>s 2 (^) = YT=i sin = m/2 has been used. 
The stresses for a disclination have already been computed l43l with result 
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The total stresses are 
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where the external pressure II is given by n = ((A - 1/2)8 + jj T m)G. Using 
Eq. 03 for the modified disclination charge s e ff = s — j^m, one finds 



F = S -^lK R 2 + ^R 2 
32vr U 8irB 



(53) 



Therefore, perfect screening implies s e ff = 0, 



b 



(54) 




m 



and there is no external pressure II = 0, the system of disclination plus 
grain boundary does diverge quadratically with the system size. It can also 
be proved that within the same approach, the term linear in system size 
vanishes also. It should be emphasized that Eq. |M]is a linear order result. 
One should expect that if all orders were included, some modifications should 
occur, among them formula Eq. [^should become Eg. I1UI 

C Energy of a Grain boundary of dislocations 

The energy of a linear grain boundary of n dislocations with Burgers vector 
perpendicular to the grain is a well known result [303 • Here we just outline 
the aspects important for this paper. The total energy is given by 



where a xy is the stress tensor of the dislocations in the xy direction. The 
final result of the energy is 



for small x, (dislocations separated a distance D >> a), the usual formula 
H(x) = — log(x) + 1 follows. If x is not small enough, the formula Eg 1571 
must be used. 




(55) 




(56) 



where a is the lattice constant and 




(57) 
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D Relations of base vectors 

If ei and e2 are the triangular denning the triangular lattice the following 
relations follow. 

£4 & < b = 3^\ (58) 

b 

where the sum is over all nearest neighbors b. It also follows 

E <b<b4 b e l ab = l(S ij S kl + FS* + W) (59) 
b 

and 

E ei b eib4b + Ab + i = l(S ij 5 kl + R U R kj + R lj R ki ) , (60) 
b 

where = cos(7r/3)5 1 - 7 + sin(7r/3)e* J . Similarly, 

E = j^tf * + B»B? k + W) , (61) 

b 
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